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Abstract 

We continue the study of the tensor-four-scalars theory which is a 
modification of general relativity. We include normal matter by ap- 
plying the displace, cat, and reflect method to our previous vacuum 
solutions with dark halo. The resulting disklike solution has physical 
and unphysical features: the matter density for large radius is propor- 
tional to the forth power of the circular velocity in agreement with the 
Tully-Fisher relation; but the radial pressure is negative. 
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1 Introduction 



In [1] and [2] we have studied vacuum solutions of the tensor-four-scalars 
theory. These solutions have the interesting property that the behavior of 
the corresponding circular velocity V{r) (rotation curve) can be arbitrarily 
given. As a consequence we find solutions with asymptotically constant 
velocity for large r [2] as observed in the dark halo of galaxies. 

The next step is the inclusion of normal matter in the theory. There 
exists an elegant method to construct solutions with a material disk from 
vacuum solutions. This is the so-called displace, cut, and reflect method 
which goes back to Kuzmin [3] and since then was used and modified by 
many authors (see [4] and references given there). The method is most 
simple if the matter is located in a two-dimensional plane which simulates 
a spiral galaxy. Then, if the plane is z = in cylindrical coordinates, the 
matter density is singular ~ 5{z) and the metric g^u has jumps in its first 
derivatives. The mathematical basis in this situation is Taub's theory of 
distribution valued curvature tensors [5]. We review the essential parts of 
this theory in the next section and give short derivations of the relations 
which we need to construct our disk solutions with dark halos. 

In section 3 we apply the displace, cut, and reflect method to the vacuum 
solution with dark halo which was derived in [2]. We obtain the surface 
energy-momentum tensor expressed by the metric functions a(r) , b(r) , c(r) . 
In section 4 we discuss this tensor in the interesting case of flat rotation 
curves which one attributes to a dark halo. The most interesting result is 
that the energy or mass density in the flat region is proportional to the 
forth power of the flat circular velocity Va at (R). We argue that this is in 
accordance with the Tully-Fisher relation. [8] [9]. On the other hand the 
radial pressure comes out to be negative which is unphysical, while the 
azimuthal pressure is correctly positive. 

2 Taub's theory of distribution valued curvature 
tensor 

We specialize Taub's general setting [4] to our situation of a planar disk. 
Let S be a three-dimensional surface in 4-space where the metric tensor 
is continuous but has finite jumps in the normal derivatives; the derivatives 
in the tangential directions are assumed to be continuous. In an admissible 
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coordinate system let S be described by the equation 

<p{x) = (2.1) 

and have the normal vector 

», - 1- 

Then the finite discontinuities in the first partial derivatives of g^ v are given 

W^L-^L—V, (2.3) 

where + and — mean the limiting values from both sides of S. This follows 
from the decomposition of the gradient into normal and tangential compo- 
nents. The corresponding jumps in the Christoffel symbols then are 

2[T%] = npb« + n 1 b c l-n a b M . (2.4) 

The Ricci tensor 

Rfiu = d a v® u — d u r® a + r^r^ — r^r a/i (2.5) 

contains derivatives of T, consequently the finite jumps lead to singular 
contributions proportional to the delta distribution 5s with support on S 
according to 

= [T« u ]np5 s . (2.6) 

Then it follows from (2.4) that 

-Raising = ^{-n a n a b^ + rf-b^Uv + n a b av n^)bs, (2.7) 

with 

b a p = b a p-\b8%, b = g^ u . (2.8) 

This is in agreement with eq.(2.14) of Taub, note that his convention for the 
Ricci tensor is the negative of ours (2.5). 

In the field equations these singular distribution must be compensated 
by a distribution valued energy-momentum tensor 



1 
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where 



(Rfiv - Rising = KtpvSs, (2.9) 



R = <f*R a p, « = ^. (2.10) 
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If the jumps of the normal derivatives of g^ u are known, t^ u can be 
calculated from (2.7) and (2.9): 

-(^-^)(^-^))b,r, (2.H) 

where n 2 = n a n a . This agrees with eq.(6-2) of Taub. This singular contri- 
bution (2.11) must be added to the regular energy-momentum tensor which 
renders the field equations fulfilled outside of the surface S. 



3 Thin disk with dark halo 



We study a simple model of a spiral galaxy by assuming that all normal 
matter is concentrated in the plane z = with a singular density ~ 5 1 (z). 
Outside this plane which is our surface S of Sect. 2 we have "vacuum" with 
a dark halo as it is described by the spherically symmetric solution of the 
tensor-four-scalars theory [2]. The corresponding line element in spherical 
coordinates (t, r, <I>) reads 



d s z = e a dt 2 - e b dr 2 - r 2 e c {d$ 2 + sin 2 $d<$> 2 ). 



(3.1) 



The metric functions a(r),b(r), c(r) are given by the circular velocity squared 
V 2 (r) = u(r) according to 



K n 



1 + J2U 



1 + 7i u 
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7i it 



1 + 72^ 



1 1 + 2u + 7U 2 



r 2 u 2 



1 + 7i u 



1 + 72^ 



with 



71 



7 



72 



(3.2) 

(3.3) 
(3.4) 

(3.5) 



1 + VT=T 

Here 7 < 1 is a parameter which measures the deviation from general rela- 
tivity and K a ,K c are positive constants of integration. 
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To have a simple representation of the plane z = and the corresponding 
delta-measure we go over to cylindrical coordinates (t, R, z, $) 

R 

r 2 = R 2 + z 2 , z = rcos$, sintf = — . (3.6) 

r 

Then the metric (3.1) assumes the following non-diagonal form 

ds 2 = g fiu dx' 1 dx 1 ' 

with 

a R h % r R r Z h 

9oo = e , 5n = s-e ~e , 522 = re ^e" 

512 = 521 = -2^(e b - e c ), 533 = -flV. (3.7) 

For simplicity we still write r, but our admissible coordinates are x 1 = 
R,x 2 = z. We also need the inverse 

5°° = e- 5 U = f, <7 22 = ^ 

5 12 = "^=5 21 , 5 33 = -, (3.8) 

-O 533 

where the determinant D is equal to 

.R?z 2 
r 



D = 511522 - (512) 2 = e b+c - 3—(e b - e c ) 2 . (3.9) 



To construct the metric with the material disk we apply the widely used 
displace, cut, and reflect method. Following the procedure of Voigt and 
Letelier [4] we take the metric (3.7) in the half space z > d > 0, displace it 
to z = and reflect it for z < 0. This produces the finite jumps in the z- 
derivatives of g^ u . The whole procedure is equivalent to the transformation 
z — > \z\ + d and setting z = at the end.. 

For the calculation of the energy-momentum tensor from (2.11) we 
need the normal vector = (0, 0, 1, 0) = 6 2 and 

n t/ = <f' i n li = <f 2 , n»n v = g 22 . (3.10) 

The jumps (2.3) in the normal derivatives are equal to 

, r 1 r%)0 dri r , zi , 2d 
boo = [goo,2\ = [-^rg^\ = [9oo-\=9oo — 
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bn = [ 9 n, 2 ] = 9^ - ^ (3.11) 
j r , , 2d 

012 = 5l2,2j = 9l2~ 

r 

, 2d 4d 6 

022 = 522,2 =922 e 

r r 

033 = [933,2 J = fey > 

where the prime always means the partial derivative with respect to r keep- 
ing z and R constant. Now from (2.11) we obtain 

2^oo = -9 22 (ooo - 9oo(b - ^^W) (3.12) 

with 

6 = gTb^ = g°Xo + 9 U bii + 2g 12 b 12 + g 22 b 2 2 + g 33 b 33 . 
After some cancellation of terms we have 

2Ktoo = goo(Db 11 + ^-b 33 ). (3.13) 
In the same way we get 

2«tn = g 22 gn(g°Xo + g 3 % 3 ) (3.14) 
2«ti2 = g 22 gi2(g°Xo + g 3 % 3 ) (3.15) 
2^22 = g 22 (^)\gn((g°Xo + g 33 b 33 ) + hi{D~* - 1)(1 - g 12 )) (3.16) 

2k<33 = 9339H (g 22 g°\o + D^bn). (3.17) 
We want to calculate the energy density 

hi™" + ■£:*»)■ (318) 



Using 



we obtain 



2k V Dg 33 



\ Or r £ 
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In the same way we obtain 

t{ = -±^e-c(a' + c') (3.20) 

*§ = -ii^- 08 ^ + R2 *- b - Cd r& + (3-21) 

Here we have to put z = everywhere which gives r 2 = R 2 + d 2 . As a 
consequence, 512 vanishes so that t\2 and £22 also vanish. 



4 Discussion of the results 

Now we must specify the circular velocity squared u(r) in order to fix the 
metric. We are particularly interested in the case of an asymptotically flat 
circular velocity which in the usual terminology corresponds to a "dark 
halo" . Therefore we assume u(r) of the form 

u(r) = u^ t + ^ + 0(r- 2 ) (4.1) 

for large r. Then it follows from (3.2-4) 

e a = K a + 0{r- 1 ), e b = ^+0(r" 5 ) (4.2) 



r 

L 

r 

where by (3.4) 

,-2 _ TA-4 



e c = 5 + 0(r- 3 ) (4.3) 



Lc ~ «££ = ^f at 4 . (4.4) 
Using this in (3.19) the leading order comes from the last term 

*0 = ^^(1 + O(i?- 1 ))- (4.5) 

This is proportional to the density of normal matter because we consider a 
static energy-momentum tensor. Taking (4.4) into account we find that 

*0~4at~<at(«) (4-6) 

for large R. This is in accordance with the baryonic Tully-Fisher relation for 
galaxies [8] [9] , which states that the total baryonic mass M is proportional 
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to Vg at . In fact, we will soon see that the contribution of the inner part 
R < Ri of the disk can be made arbitrarily small compared to the outer 
part between R\ < R < R2, say. We emphasize that M is obtained from 
t|j by integrating with the euclidean surface measure RdRdQ because this 
is what astronomers are doing when they determine M from luminosity 
measurements. 

For small R the circular velocity behaves linearly V(R) ~ R so that 

u(r) = u 2 r 2 + 0(r 3 ). (4.7) 
This gives the following behavior of the metric functions (3.2-4) 

e a = K a + 0(r 2 ) (4.8) 

e fe = ^ + 0(r" 5 ), e c = ^ + 0(r- 5 ). (4.9) 

Note that everything remains finite for R — >■ because we have r 2 = R 2 + 
d 2 > d 2 . Now the leading order in (3.19) comes from the first term which 
gives the positive contribution 

+ (4.10) 

This can be made arbitrarily small by choosing R\ small or d big enough. 

The simple displace, cut, and reflect method often gives unphysical re- 
sults for the radial pressure t\ [4]. Our solution (3.20) suffers from the same 
defect. The leading order in the outer region R\ < R < R2 corresponding 
to (4.1) comes out to be positive 

Since this is the negative radial tension —Pr the latter is negative which is 
impossible for normal matter. On the other hand the azimuthal pressure t| 
has the correct sign 



which gives a positive azimuthal tension 

Summing up, the simple thin disk model of a spiral galaxy with dark 
halo gives interesting results as far as the dark matter problem is concerned. 
For a more realistic description a better model is needed. All our conclusions 
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remain true if we consider the classical case of general relativity, i.e. 7 = 0. 
So this good old theory in the non-standard gauge [2] still remains an option 
to understand the dark matter problem. But we have some indication that 
the future better theory only works with the tensor-four-scalars theory 7 > 
0. 
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